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Abstract. — Taking advantage of the advances in array detector technology, an imaging polarimeter (IMPOL) 
has been constructed for measuring linear polarization in the wavelength band from 400-800 nm. It makes use of a 
Wollaston prism as the analyser to measure simultaneously the two orthogonal polarization components that define 
a Stoke’s parameter. An achromatic half-wave plate is used to rotate the plane of polarization with respect to the 
axis of the analyser so that the second Stoke’s parameter also can be determined. With a field of view correponding 
to about 30 X 30 mm^ for a (2)1.2 m, f/13 telescope, a sensitive, liquid-N 2 cooled CCD camera as the detector and a 
built-in acquisition and guidance unit, the instrument can be used for studying stellar helds or extended objects with 
an angular resolution of ~2" . The instrumental polarization is less than 0.05% and the accuracies of measurement 
are primarily limited by photon noise for typical observations. 
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1. Introduction 

The advances in two dimensional array detector technol¬ 
ogy in the optical and near infrared wavelength bands 
have made new kinds of imaging astronomical observa¬ 
tions feasible. Astronomical polarimetry is one field which 
has gained tremendously from these developments. The 
limitations of using aperture photometry for polarimetry 
were so severe that any serious study was rendered time 
consuming and difficult. On the other hand, imaging po¬ 
larimetry with its capabilities for multiplexing, simultane¬ 
ous sky measurement, seeing-limited resolution etc. offer 
great advantages over aperture polarimetry. Astronomers 
have recognized this potential and have developed several 
new observation techniques in the optical (eg. Scarrott 
1991; Jannuzi et. al. 1993; Jarrett et. al. 1994; Wolsten- 
croft et. al. 1995; Simmons et. al. 1995) and near infrared 
(eg. Kastner & Weintraub 1994; Moore & Yamashita 1995; 
Weintraub et. al. 1995) wavelenghts to study phenomena 
in a variety of Galactic and extragalactic astrophysical 
objects. In this paper, we report the design and construc¬ 
tion of an imaging polarimeter (IMPOL) which uses a 
cooled CCD array as detector. It was developed at the 
Inter-University Centre for Astronomy and Astrophysics 
(lUCAA), INDIA. The principle of the instrument (Sen 
& Tandon 1994) is based on a combination of ideas sug- 
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gested by Ohman (1939) and Appenzeller (1967). An in¬ 
strument of this type has been constructed at the Univer¬ 
sity of Durham and has been in use for some time now 
(Scarrott et. al. 1983). Section 2 is a description of the 
instrument - the different subsections dealing with vari¬ 
ous aspects ranging from design guidelines to instrument 
control and user-interface. The dominant sources of errors 
in the measurement are investigated in Sect. 3, while an 
estimate of the performance of the instrument under two 
typical observing conditions is given in Sect. 4; Section 4 
also contains an estimate of the performance of the acqui¬ 
sition & guidance unit. In Sect. 5 we discuss the observa¬ 
tional procedure and the data-analysis method. Section 6 
contains results of the commissioning tests of the instru¬ 
ment. The last section (Sect. 7) contains the concluding 
remarks. 

2. The Instrument 

The following set of scientihc and technical guidelines were 
adopted as basis for the design of the instrument, (i) It 
should be capable of observing faint extended objects like 
reflection nebulae, accretion disks, dusty active galaxies 
etc., with accuracy limited by photon noise and resolu¬ 
tion limited by seeing. This demands a typical field of 
view of a few arcminutes, an optics which is well-matched 
with the telescope to minimize loss of light and a detec- 
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tor of sufficient sensitivity. Also, the effects of atmospheric 
scintillation should be eliminated and instrumental polar¬ 
ization minimized so that such measurements of low level 
polarization could be made, (ii) Since the objects of inter¬ 
est are very often quite faint, the instrument should have 
an acquisition and guidance (A&G) system which allows 
pointing the telescope with an accuracy of a few arcsec- 
onds and tracking to better than 1" over long periods, 
(iii) Multi-wavelength observations should be possible in 
various wavelength bands in the optical and near-IR re¬ 
gions. (iv) The entire instrument consisitng of the optics, 
A& G unit, associated electronics etc. should be a self- 
contained unit which can be easily mounted on the tele¬ 
scope, with only electrically isolated communication links 
to the computers for instrument control and data acquisi¬ 
tion. (v) The cost of the instrument should be minimized 
by using standard optical and electronic components and 
by avoiding over-specification as far as possible. 


2.1. Principle of the instrument 

Figure 1 is a schematic representaion of the optical ar¬ 
rangement of IMPOL. The basic idea behind this arrange¬ 
ment is to use a Wollaston prism with its axis normal to 
the optical axis of the system, as the analyser to convert 
the linear polarization in the incoming light into relative 
intensity of two orthogonally polarized beams (the ordi¬ 
nary and the extraordinary), separated by a small angle 
of 0.5°. This measurement is sufficient to define one of the 
Stoke’s parameters Q or U. The other Stoke’s parameter 
can be determined by rotating the plane of polarization 
of the incoming light relative to the analyser through a 
known angle and taking another measurement. This is 
accomplished by introducing a half-wave plate, with its 
fast-axis normal to the optical axis of the system, be¬ 
fore the Wollaston prism. When the half-wave plate is ro¬ 
tated through an angle a, the plane of polarization rotates 
through an angle 2a. At this new position of the half-wave 
plate another measurement on the orthogonally polarized 
beams can be made to determine the second Stoke’s pa¬ 
rameter as well. It is easily seen that, for this arrangement 
the intensities of the extraordinary and ordinary beams 
(/e and lo) are given by 



Fig. 1 Schematic of the IMPOL optical layout 


Ie{a) = ^-b JpCOS^(6»-2a) (1) 

loia) = ^-b/pSin^(6» - 2a) , (2) 

where /up & Ip are the intensities in unpolarized and po¬ 
larized condition respectively in the incoming beam; 6 a 
are the position angles of the polarization vector and the 
half-wave plate fast-axis respectively with reference to the 
axis of the Wollaston prism. Since the angle 9 is conven¬ 
tionally measured with respect to the celestial north-south 
axis (and increasing counter-clockwise), the axis of the 
Wollaston prism is kept aligned to it. We define the ratio 

¥ -1 

R{a) = -f -= pcos{29 — 4a) , (3) 

To + ^ 

where p = r \r , is the fraction of the total light in fin- 

-‘up 

early polarized condition. This ratio reduces to the nor¬ 
malized Stoke’s parameters y and y for a = 0° and 
a = 22.5°. In practice, additional measurements are made 
at two more values of a, namely 45° and 67.5°, for reasons 
explained in Sect. 5. 

The half-wave plate and the Wollaston prism are 
placed in between a field lens-camera lens combination, 
which reimages the telescope focal plane on to the main 
GGD with a reduction factor of about 2.7; the field lens 
reimages the telescope aperture on the half-wave plate and 
the fight reaches the camera lens without any vignetting. 
As shown in Fig. 1 each point in the telescope focal plane 
produces two images on the GGD, corresponding to the 
ordinary and the extraordinary beams. In order to avoid 
overlap of the images of adjacent points, for observations 
of extended objects, a grid of parallel obscuring strips is 
placed at the focal plane of the telescope. The width of 
and spacing between the strips are chosen in such a way 
as to avoid the overlap of the ordinary and the extraor¬ 
dinary images on the CCD. Four 0.3 mm diameter holes 
are provided at the four corners of the grid and are used 
to focus the surface of the grid, which coincides with the 
focal plane of the telescope, on to the surface of the CCD. 
The grid is made of black dielectric material to avoid po¬ 
larization of the stray light arising due to reflections from 
its edges. Besides, the edges are made slanted (Fig. 1) to 
prevent vignetting of the telscope beam. 

2.2. Acquisition and Guidance Unit 

A schematic illustration of the Acquisition and Guidance 
(A&G) unit is shown in Fig. 2. The right-angle prism, the 
field lens (LI) and a relay lens (L2) are held in position 
inside a probe which in turn is mounted on a linear XY 
translator stage. This arrangement allows the probe to be 
positioned anywhere within a 40 x 100 mm^ rectangular 
area (corresponding to about 45 square minutes of arc for 
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Fig. 2 Schematic of the Acquisition & Guidance unit 


a ©1.2 m, f/13 telescope) close to the main field. The lens 
L2 reimages a ©12 mm area of the telescope focal plane on 
to the tip of an image quality optical fibre which, in turn, 
conveys it to the photocathode of an image intensifier. 
Two Nikor lenses (L3 and L4) focus the intensiher output 
on the surface of the CCD of the A&G unit. 

For pointing the telescope, atleast one star brighter 
than about 15*^ mag. in visual and close to the main field 
is identified and its right ascension and declination values 
are entered in to the control computer along with those 
of the main object. Once the translation stages are moved 
to the requisite offset coordinates by the computer, short 
exposures are taken with the CCD of the A&G unit to 
adjust the telescope pointing such that the guide star ap¬ 
pears at the center of the CCD. Now the telescope has 
been pointed and auto-guiding is started. During auto- 
guiding the control computer takes successive exposures 
with the CCD and hnds the centroid of the guide star 
image. Error signals are then generated accordingly and 
transmitted to the control system of the telescope. 

Table 1 is a list of the optical components used in IM- 
POL along with their sources and part numbers. The field 
lens, half-wave plate, Wollaston prism and the filters were 
unmounted when purchased. 


2.3. Instrument control, Data-acquisition and User- 
interface 

A mixture of digital and analogue electronics hardware 
and a layered architecture software have been used to per¬ 
form the instrument control and data-acquisition opera¬ 
tions. Apart from the control computer (PC 486), all the 
electronics is mounted on the telescope along with the in¬ 
strument. The upper three layers of the software reside 
in the control computer, while the lower two layers are 
downloaded into the local memory of the instrument at 
the time of start-up and initialization. 



Fig. 3 Various blocks of the IMPOL control system. The 
block marked “T” represents the XY stages for position¬ 
ing the probes. 


2.3.1. Electronics 

The electronics assembly of the instrument consists of 
three parts - (i) the positioning systems for the half-wave 
plate and the A&G unit probe, (ii) the CCD camera for 
the A&G unit and the associated electronics and (hi) the 
main CCD camera with its own control electronics and 
host computer (PC 486) for exposure control and data- 
acquisition. A block schematic of the IMPOL control sys¬ 
tem is shown in Pig. 3. The half-wave plate and the A&G 
unit probe are positioned using stepper motors which are 
controlled with signals generated by a 8051 microcon¬ 
troller based card. These digital control signals are fed to 
translator circiuts which convert them to power-amplified 
phase signals required to drive the stepper motors. Sen¬ 
sors are used on the half-wave plate mount and the A&G 
unit XY positioning stages for fixing the zero posotions. 
The controller card, in turn, receives “command” packets 
from the computer for instrument control. 

The detector for the A&G unit is the commercially 
available ST-6 CCD camera of the Santa Barbara Instru¬ 
ment Group (SBIG). This camera has a controller unit 
of its own which executes specific tasks under the com¬ 
mand of the control computer. (The SBIG have, in a pri¬ 
vate communication, graciously provided us the instruc¬ 
tion set for their controller which we have used in our 
custom-software for instrument control.) The “command- 
response” exchange between the control computer and the 
instrument takes place over an optically isolated asyn¬ 
chronous serial communication link (RS 422 standard). 
The ST-6 controller has a second serial communication 
port (AUX port) so that it can be used as a gateway for 
communicating with other peripherals which use the same 
type of communication link. The 8051-based controller 
card, described above, is connected to this AUX port and 
the communication between it and the control computer 
is routed through the ST-6 controller (see Fig. 3). 
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Table 1. Optical components used in IMPOL 


Component 

Source 

Part number 

Description 

Field lens 

Karl Lambrecht 

322305 

(2)50 mm, eff. EL 300 mm, 
AR coated, achromat 

Half-wave plate 

Karl Lambrecht 

WPAC 2-22-BB400/700 

achromatic, 22 mm clear aperture, 
broad band single layer AR coated 

Wollaston prism 

Bernard Halle 

PWQ 30.30 

quartz, 29x29 mm clear aperture, 
multi layer AR coated 

Filters 

Andover 


B, V and R, 040 mm, 5 mm thick, 
AR coated, image quality 

Camera lens 

Nikon 


eff. EL 85 mm, f/1.8, A.F. 

Right-angle prism 

Melles Griot 


12x12 mm 

Lenses LI 

L2 

Edmund Scientific 


eff. EL 25 mm, f/2.0 
eff. EL 6 mm, f/2.0 

Optical fibre 

Dolan Jenner 

ISO 636 

6x6 mm, coherent glass fibre, 
36 inch long, N.A. 0.66, 

Image intensifier 

Philips 

XXI500 

014 mm, resolution 36 Ip/mm 

Lenses L3 & L4 

Nikon 


eff. EL 50 mm, f/1.2 


The main CCD camera for data-acquisition is built 
around a liquid N 2 cooled, 385x578 pixel, front-side il¬ 
luminated, phosphor coated CCD chip from EEV. The 
operations of the CCD camera like shutter control, bin¬ 
ning mode selection etc. are carried out by another 8051- 
based controller card which communicates with the data- 
acquisition computer over a second asynchronous serial 
communication link. Once the exposure is over the image 
pixal data are transmitted to the computer over a high 
speed cable using dedicated circuits. The important pa¬ 
rameters of this CCD camera, which was also developed 
at lUCAA, are listed in Table 2 (for more details refer to 
Deshpande and Cadre 1994). 

2.3.2. Software 

The instrument control software is made of five layers, of 
which none but the top most one is visible to the user. 
This layer provides the user-interface, reads the system 
configuration files and initializes the system. It is at this 
level where the sequence of operations is decided, depend¬ 
ing upon the instructions given by the user and corre¬ 
sponding service calls to the immediate lower level are 
generated. The second layer provides these services mostly 
by creating a “command-l-parameter(s)” packet depend¬ 
ing on the service requested and the current system con¬ 
figuration and then passing this packet on to the third 
software layer. The communication protocol between the 
control computer and the microcontrollers at the instru¬ 
ment is defined in the third layer which sents the packet 
to the ST-6 controller or through it to the 8051-based 


Table 2. The parameters of lUCAA CCD Camera 


Parameters 

Value 

CCD make 

EEV CCD02-06 series 

CCD chip size 

385(H) x578(V) 

Pixel size 

22 X 22 /im 

Active area 

12.7 mm x 8.5 mm 

No. of amplifiers 

I 

Quantum efficiency 

blue : 20% 
yellow : 50% 
red : 60% 

Readout speed 

32 fj,s per pixel 

Acquisition & Display time 

12 s for full frame 

Read noise (total) 

8 e“ rms 

Gain 

5 e- / ADU 


controller for positioning operations. Upon receiving a 
“command-hparameter” packet, either of these controllers 
parses it and checks for errors. If the command is valid 
and there are no errors it immediately sends back an ac¬ 
knowledgement signal. Otherwise it sends back a corre¬ 
sponding error signal. The result of this communication 
process is conveyed back to the highest layer for appropri¬ 
ate follow up. It is to be noted that at this stage the top 
most layer does not have information whether its service 
request has actually produced the required effect or not. 
To obtain this information it has to generate another ser- 
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vice call requesting the current status of the appropriate 
sytem parameters. 

The 8051 microcontroller has a built-in interpreter for 
BASIC. The processes of packet reception, parsing, error 
checking and acknowledging are carried out by a BASIC 
program which resides in the memory of the microcon¬ 
troller forming the fourth software layer. A set of assem¬ 
bler routines, each doing a specific operation, forms the 
fifth and lowermost software layer. The BASIC program 
calls one of these assembler routines to execute the par¬ 
ticular task requested by the command. 

This system of control encashes on the fact that while 
each individual operation is to be carried out fast, the 
interval between successive operations are comparable to 
typical user response time scales. Though the layered ar¬ 
chitecture appears complicated, it is in fact extremely easy 
to implement and debug and provides great flexibility to 
alterations. 

3. Sources of errors 

In addition to the fundamental limit due to photon 
statisitcs, the measurement of polarization is affected by 
several factors beginning from the atmosphere down to the 
detector and the data-analysis procedure. In the following 
four subsections we look at, in some detail, the important 
sources of errors and the techniques used to minimize their 
effects. 

3.1. Atmospheric Effects 

Since polarization measurement involves taking multiple 
exposures, variations in atmosperic transparancy and scin¬ 
tillations can affect the measurement process. In IMPOL, 
two orthogonal polarization components are measured si¬ 
multaneously and the Stoke’s parameter is obtained from 
the ratio of the fluxes in these components (Eq. (2)). Since 
the atmosphere is not birefringent, this eliminates the ef¬ 
fects due to atmospheric scintillation, or for that matter 
any effect which changes both the polarization compo¬ 
nents by the same factor-like variations in the effective 
exposure times of observations, presence of thin clouds 
etc. 

3.2. Distortions in polarization 

As the measurements depend on detecting changes in the 
polarization vector caused by rotating the half-wave plate, 
any small instrumental polarization occuring after the 
half-wave plate can be easily eliminated. All the optical 
elements are anti-reflection coated to minimize the po¬ 
larization efffects due to reflection at their surfaces, and 
care is taken to minimize the stray light, reflected from 
the walls etc., reaching the detector. In particular, all the 
aperture stops in the beam path are chosen to be non- 
metallic, including the grid in the focal plane, so as to 


avoid polarization of light scattered from these (Posperge- 
lis, 1965). The mounts for field lens, the half-wave plate 
and the Wollaston prism have been designed to minimize 
stress-birefringence due to differential thermal expansion. 

The retardance introduced by the half-wave plate 
could deviate from 180° either because of the finite an¬ 
gle of incidence or because of the chromatic effects. If 
the beam is incident at a small angle i, the maximum 
change in retardance is given by (derived from expression 
in Serkowski 1975) 


CTr — TT X 


4no 



( 4 ) 


where Uq and Ue are the refractive indices of the mate¬ 
rial for ordinary and extraordinary rays. For the aperture 
used, the maximum angle of incidence is about 5° so that 
(Tt ^ 0.01 rad. The chromatic effects give ctt < 0.18 rad 
for the wideband and Ct < 0.14 rad for the V-band. It 
can be shown (Serkowski 1974) that the depolarization Up 
due to an uncertainty of Or rad in the retardance is to 
the lowest order given by dp ~ pa^, p being the fractional 
polarization. Therefore even for the wideband dp ^ 0.03p. 
Further, circular polarization in the incident light is con¬ 
verted to a linear polarization of magnitude Vdr, where V 
is the circularity parameter. For wideband observations, 
conversion to linear polarization would be about 0.18F, 
but in typical observations this does not pose a serious 
problem because the circular polarization in the incident 
light is usually less than the linear polarization. 

Another chromatic effect is the change in position an¬ 
gle of the half-wave plate fast axis with wavelength. Al¬ 
though this does not produce any error in the measure¬ 
ment of fractional polarization p, it does render the mea¬ 
surement of position angle 9 erroneous. A solution to this 
problem is to observe standard stars using narrow band 
filters and produce a smoothed curve which gives the re¬ 
quired position angle correction for each wavelength. The 
achromatic half-wave plate used in IMPOL does not pro¬ 
duce any appreciable dispersion in its fast-axis position 
angle over the wavelength range of interest. 

The uncertainty involved in the positioning of the half¬ 
wave plate leads to an error in the measurement of po¬ 
larization. From Eq. (2), we can see that a positioning 
uncertainty of da produces a maximum error in the mea¬ 
surement of linear polarization given by dp ~ pda- Thus 
for an error of 0.1°, of the postioning system, this results 
in dp cs 0.002p. 


3.3. Flat held errors 

The half-wave plate is the only moving component in the 
optical train and hence it is relatively easy to make sure 
that the images do not shift on the detector. Further, as 
the analyser is fixed with respect to the detector, the ori¬ 
entation of the ordinary and extraordinary polarizations 
also remain fixed with respect to it, independent of the 
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polarization vector of the incident beam. This renders it 
easy to estimate the flat field correction factor from the 
set of observations at the four positions of the half-wave 
plate. These points are further explained in Sec. 5. 

3.4. Photon noise 

From the above discussion we find that all the sources of 
errors described do not contribute more than 0.05p to the 
error in the measurement of p. In the next section, we es¬ 
timate the errors under typical observing conditions that 
result from photon noise alone. Comparing the two, we 
conclude that the dominant source of error in the typical 
measurement of polarization is photon noise. From equa¬ 
tion (2) it can be shown that the variance due to photon 
noise in the measurement of each of the Stoke’s parame¬ 
ters R{a) is given by 

where Iq and R are the number of photoelectrons and k is 
the ratio of the flux from the background to that from the 
source. For small values of p, /q and /e are approximately 
equal and denoting I = Io + le, the standard deviation of 
p can be written as 


distributed and might have bias (Clarke et. al. 1983). Be¬ 
sides, for photon-noise dominated measurements, the pos¬ 
itive kurtosis of the distribution will lead to erroneous es¬ 
timates of the confidence levels unless more than a few 
thousand photoelectrons are collected. Thus, in order to 
arrive at an optimum procedure, it is essential to carefully 
study the nature of the dominant sources of noise in the 
measurement and their effects on the analysis (see Sec. 5). 

4. Performance estimates 

An estimate of the performance of the instrument is 
made in the following two sub-sections. In the first part 
we calculate the accuracy of polarization measurement 
as limited by photon noise for two typical types of ob¬ 
servations. The second part deals with the performance 
of the A&G unit. In both cases, it is assumed that 
the observations are made in the V-band (bandwidth 
AA ~ 100 nm, Vq = 1000 photons per cm^ per s per 
A) using a (l)l.2 m (area A ~ 10^ cm^) telescope having 
an f/13 beam, that the brightness of the sky is 20 mag. 
per sq. arcsec and that the software aperture used for pho¬ 
tometry has about 30 pixels, ie. about 30 sq. arcsec. 

4.1. Polarimetry 


ap ~ 100 X , 

where Ib = kl. Similarly, we can also show that 

CT p 

ag ~ 0.5 X — rad . 

P 


( 6 ) 

(7) 


It is worth noting here that p, as defined above, is a 
positive definite quantity and follows the Rice probability 
distribution given by 


F{p,Po) 


■7o(^) 


( 8 ) 


Here, po is the true value of fractional polarization being 
estimated by p and Iq is the modified Bessel function of 
order zero. Since p is a biased estimator, several schemes 
have been suggested (Simmons & Stewart 1984 and ref¬ 
erences therein) for debiasing, but none of them is fully 
satisfactory. For values of larger than about 4, almost 
all the debiasing schemes agree and reduce approximately 
to the relation _ 

Po = 

where po is the “debiased” estimate of po- Thus it is de¬ 
sirable to work with the normalized Stoke’s parameters 
as far as possible and use the quantities p and 9 only to 
present the final results. However, in general, the normal¬ 
ized Stoke’s parameters themselves might not be normally 


In the following discussion the combined effect of the at¬ 
mosphere, the telescope, the optics and quantum efficiency 
of the detector etc., is found to result in an effective trans¬ 
mission of ?7 ~ 3.2%. 

A classical example of polarimetric observations in¬ 
volves study of dark molecular clouds with the light of 
stars shining behind their periphery (eg. Elvius 1970; Joshi 
et al. 1985; Kane et al. 1995). Such studies give valuable 
information regarding the properties of dust in the cloud 
and the magnetic field in its vicinity. Typically the po¬ 
larization of these sources is in the range 1-5 %. Assum¬ 
ing stars of apparant visual magnitude TOv = 15, we see 
that the total number of photo-electrons collected, due to 
the star (J) and the background (Ib), within the software 
aperture, for an exposure of 15 minutes are 2.9 x 10® and 
8.6 X 10"^ respectively. Therefore, the corresponding error 
in polarization measurement is given by 


Op = 100 X 


V2.9 X 10® -k 8.6 X 104 
2.9 X 10® 


0 . 21 % . 


( 10 ) 


If the observations are made with a filter which covers 
both the V and R bands (AA ~ 200 nm), there will be 
an improvement in the signal-to-noise ratio by a factor of 
72, giving cTp ~ 0.15%. 

Study of extended objects like reflection nebulae, is 
another field in which polarimetric observations are use¬ 
ful. Polarization in the light scattered from these clouds 
can range from a few percent to as much 30%. Assuming 
that the observation involves a reflection nebula of the 



Page Header Title 


7 


same surface brightness as the background sky, the inte¬ 
grated photoelectron count from each, for an exposure of 
15 minutes, will be 8.6 x 10^ and hence 



12 14 16 18 20 

visual magnitude 


Fig. 4 The errors due to photon noise alone, as calculated 
in Sec. 4.1, for wideband measurements of fractional po¬ 
larization p are shown as a function of the brightness of 
the source. The solid curve refers to the case of stars while 
the dashed curve is for an extended object, for which the 
x-axis represents the surface brightness of the source in 
magnitudes per sq. arcsecs. The background is assumed 
to be 20*^ magnitude per sq. arcsec and a software aper¬ 
ture of 30 sq. arcsec has been used during data reduc¬ 
tion. The crosses represent the errors in actual polarimet- 
ric measurements made with instrument on a field at the 
periphery of the dark cloud B133. 

for a solid angle defined by the software aperture. The 
wideband observations, in this case will give an error of 
about 0.34%. 

In Fig. 4 the solid and dashed curves show how the 
error in the measurement of fractional polarization, based 
on photon statistics alone, changes as a function of source 
brightness, for the two cases discussed above, namely stel¬ 
lar sources and reflection nebulae respectively. 

4.2. Acquisition & Guidance (A&G) unit 


area has been chosen on the basis of the requirement that 
three stars brighter than ruv = 15 should be available with 
95% confidence level towards the poles of the Galaxy. 

In this case, apart from the atmospheric and telescope 
transmission losses, account has to be made for losses in 
the probe optics, optical Hbre and the quantum efficiency 
of the intensiher. Together this amounts to an effective 
transmission of ~0.1. The gain of the intensiher can be 
varied from 2000 to 50000 and the effective transmission 
of the lenses L3 and L4 and the ST-6 CCD is about 3%. 
Therefore the signal-to-noise ratio remains essentially un¬ 
changed in the stages after the intensiher. An exposure 
of 1 s for a star of V ~15 mag., gives a signal-to-noise 
ratio of about 20 in each of the 10 pixels covered by the 
image. This is sufficient to obtain an auto-guidance accu¬ 
racy better than 0.1'' with a correction frequency of about 
a second. However, in order to avoid excessive jitter, hys- 
terisis has been introduced into the system so that the 
corrections will occur only if a minimum shift of 0.2" has 
been sensed in the stellar image centroids. 


5. Observations and Data-analysis 

We have seen in Sect. 2.1 that observations for at least two 
positions of the half-wave plate are required to determine 
the three unknown parameters namely, the total intensity 
(/), fraction of light in linearly polarized condition (p), and 
position angle of the plane of polarization (0). However, 
the situation is not so simple in reality because - (i) the 
responsivity of the system to the two orthogonal polar¬ 
ization components may not be the same, and (ii) the re¬ 
sponsivity of the CCD is a function of the position on its 
surface. Due to these effects the signals which are actually 
measured in the two images (/' and I^) are given by 

/'(a) = Ieia)Fe{x,y) and 

/'(a) = Io{a)Fo{x,y) , (12) 

where Fe{x,y) and Fo{x,y) represent the effects men¬ 
tioned above; {x, y) being the coordinates on the array 
surface. If it is ensured that the ordinary and extraordi¬ 
nary images do not move on the surface of the CCD during 
different exposures, the ratio of the factors Fq and F^ can 
be estimated as 

Fo r Jo(0°) /'(45°) /'(22.5°) /'(67.5°) . 1 

F; ^/e(45°) /'(0°) /'(67.5°) /'(22.5°)J ’ ^ ’ 

by making use of the fact that a rotation of the half-wave 
plate by 45° simply leads to an interchage of the signals 
in the ordinary and extraordinary images. Now the actual 
ratio of the fluxes in the two images may be recovered as 


/e(a) fo ,, J'(a) ^ 

Jo (a) Fe I^{a) 


(14) 


The A&G unit has a field of view of about 2' and can This ratio is substituted in Eq. (2) and a cosine curve is 
be positioned anywhere within a 40x100 mm^ area at the Htted to the four values of R{a) obtained so as to make 
focal plane of the telescope close to the main field. This the best estimates of p and 9. 
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In practice, it might also become necessary to take 
several exposures at each half-wave plate position so as 
to collect sufficient number of photons to get the required 
signal-to-noise ratio. This might also become necessary if 
there are stars in the field with very large magnitude dif¬ 
ferences. Then a set of short exposures have to be taken 
to avoid the bright stellar images from saturating while 
a set of longer exposures will be required to collect suf¬ 
ficient photons from the fainter stars. In order to ensure 
that the normalized Stoke’s parameters follow a normal 
distribution as close as possible, the basic observational 
quantities, namely the fluxes measured at each position 
of the half-wave plate, should be averaged together before 
taking their ratios to give the normalized Stoke’s param¬ 
eters. However, this method has its demerits too. Firstly, 
it leaves the non-linear fitting technique with only one 
degree of freedom, to estimate the two parameters p and 
9 from the four values of R{a). Secondly, it does not al¬ 
low the estimation of errors in the measurement of the 
Stoke’s parameters, since there are not enough degrees of 
freedom. Therefore, it appears optimum to use a fitting 
technique to estimate p and 9, from the values of R{a) 
obtained from individual exposures. 

Since a grid of parallel obscuring strips is placed at the 
telescope focal plane (Sect. 2.1), only less than half the 
field is seen by the detector during any observation. Once 
a set of exposures has been taken for all the four positions 
of the half-wave plate, the telescope orientation is slightly 
changed so as to see the other part of the field which was 
previously obscured by the grid strips. Alternatively, in 
the case of stellar fields with slowly changing (in intensity 
and polarization) background, the grid may be removed 
during observations, provided the field is not too crowded. 

The images are collected by the data-acquisition com¬ 
puter during observations and later transferred to car¬ 
tridge tapes. The data analysis is carried out on UNIX- 
based workstations. A polarimetry package has been de¬ 
veloped for this purpose within the IRAF environment 
using a mixture of standard IRAF tasks, custom-made 
CL scripts and FORTRAN routines. The images are first 
subjected to the preliminary steps of bias and dark re¬ 
moval, cosmic ray and bad pixel detection, and masking 
etc. PSF fitting tasks of the DAOPHOT package are then 
used to determine accurately the centroids of the stellar 
images. The intensity estimates are, however made us¬ 
ing aperture photometry covering a diameter greater than 
2xFWHM so as to integrate more than 90% of the signal. 
(The FWHM is kept greater than 3 pixels to minimize 
the effects of undersampling and intrapixel variation of 
quantum efficiency.) 

6. Commisioning of IMPOL 

The commisioning of the instrument involved two levels of 
tests. Firstly laboratory tests were conducted on the opti¬ 
cal bench to determine the various parameters like linear¬ 


ity, gain and readnoise of the detector, vignetting, depo¬ 
larization, instrumental polarization etc. Some necessary 
modifications were made in the optical system and the 
software in order to optimize the instrument performance. 
Measurements were made with a 3000 K filament source 
in B, V, R bands, and a wide band defined by a 3 mm 
thick KG3 glass. A 100% polarized beam showed that the 
depolarization is about 0.4% in the B and V bands, less 
than 1% in the R band, and less than 6% in the broad 
band. There was no discernible rotation of the position 
angle between these bands. For R and the wideband, the 
accuracy of measurements is limited by the performance of 
the polarizer used. In order to estimate the instrumental 
polarization, a Lyot depolarizer was introduced in front of 
the artificial star and the polarization was measured. This 
experiment showed that the instrumental polarization is 
about 0.03% in the B-band and less than 0.06% with the 
wideband filter. 

Table 3. Observations of standard stars 


Star 

P 

(Jp 

9 

0-0 

f^phot 

Po 

00 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

HD 

% 

% 

0 

O 

% 

% 

O 

43384 

3.05 

0.04 

172 

0.3 

0.03 

3.0 

170 

154445 

3.65 

0.03 

90.5 

0.2 

0.02 

3.7 

90 

102870 

0.03 

0.01 

64.4 

7.8 

0.02 

0.017 

162 

39587 

0.03 

0.03 

140 

27.5 

0.02 

0.013 

20 


0.05 

0.02 

159 

12.7 

0.02 




0.09 

0.02 

34 

6.7 

0.02 




0.05 

0.02 

167 

11.3 

0.02 




The second stage of tests where conducted when 
the instrument saw the first light in February 1996, at 
the 01.2 m Gurusikhar Infra-Red Telescope facility on 
Mt. Abu, Rajastan, operated by the Physical Reearch 
Laboratory at Ahmedabad, India. Several standard po¬ 
larized and unpolarized stars were observed at this time. 
Since the standard stars are too bright, having visual mag¬ 
nitudes in the range of 4 to 7, instrument was slightly defo- 
cussed to avoid saturation of pixels in reasonable exposure 
times. 

Table 2 contains the results of the observations in V- 
band. Column 1 of the table gives the identification num¬ 
ber of the star. Columns 2-5 contain the measured values 
and Icr errors of the fractional polarization p and the posi¬ 
tion angle 9. Column 6 gives the error based on photon 
statistics alone, derived as illustrated in Sect. 3. Cols. 7 
and 8 list the published values of p and 9. 

The first two stars in the table are polarized standards 
while the last two are unpolarized ones. Comparing the 
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measured and published results for the polarized stars we 
see that there is no indication of depolarization. Neither 
is there any rotation of the position angle beyond a cou¬ 
ple of degrees. From Col. 3 and Col. 6 for the unpolarized 
stars, we see that the error in the measured polarization 
is comparable with that expected on the basis of photon 
statistics alone. In Fig. 5, the values of g = R{0°) are plot¬ 
ted against u = R{22.5°) for all the unpolarized standards 
observed in the V-band. The extremely low correlation co¬ 
efficient between q and u, is indicative of the low instru¬ 
mental polarization. The average values of q and u, give a 
value of polarization p = 0.03%. These observations show 
that for accuracies >0.05% in the measurement of polar¬ 
ization p, the performance is still limited by photon noise 
and the instrument polarization floor has not yet reached. 
The four measurements of star HD39587 are made at dif¬ 
ferent corners of the CCD frame in order to verify the 
uniformity of response. 


The data which is required to incorporate the A&G 
system was also collected during this observation run. The 
A&G system was succesfully installed during a second ob¬ 
servation run with the same telescope in late May, 1996. 
This time, several observations of stellar Helds at the pe¬ 
riphery of some dark molecular clouds were also made. 
The analysis of this data is in progress at present and will 
be reported separately. The accuracy of the polarization 
measurement can be gauged from the results for some of 
these stars shown in Fig. 4. 


7. Conclusion 


An Imaging Polarimeter (IMPOL) has been constructed 
which uses mostly standard optical and electrical compo¬ 
nents, but through careful design has been able to achieve 


q - 0.022, Sq - 0.044; u, = -0.022 = 0.070; p(q, u) - 0.063 


0.05 


0 

s 

- 0.05 


- 0.1 


- 0.15 - 

- 0.05 0 0.05 0.1 

q(%) 

Fig. 5 Normalized Stoke’s parameter q is plotted against u 
for a number of unpolarized standard stars. The correla¬ 
tion coefficient of the points is about 0.06 and the average 
values of q and u give a value of p = 0.03%. 

photon-noise limited performance. The instrument has a 
field of view of about 6.5' for a 01.2 m, f/13 telescope, 
and the detector is a sensitive, liquid N 2 cooled GGD. 
An off-axis acquisition and guidance unit capable of using 
stars as faint as Knag ~ 15 is also built in the instrument 
so that long exposures of faint extended objects like re¬ 
flection nebulae etc. can be taken while keeping the image 
fixed on the GGD face to within one half of a pixel - this 
stability of the image is a necessary condition to achieve a 
high accuracy in relative photometry between the differ¬ 
ent frames used to estimate the polarization. Observations 
of nearby standard polarized and unpolarized stars show 
that for wideband observations, there is no discernible 
depolarization and the instrumental polarization is less 
than 0.05%. Preliminary results of wideband polarimetry 
of stellar fields, using a PSF fitting technique for deter¬ 
mining the centroids of the stellar images and aperture 
photometry with a radius of 3 pixels to derive the inten¬ 
sities of the individual images, give close to photon-noise 
limited accuracy of 0.15% for Vmag = 15 stars with about 
an hour of total exposure time. For extended objects with 
brightness about 20th mag. per square arcseconds and a 
background of about the same brightness, it is estimated 
that an accuracy of about 1% is possible with 60 minutes 
of total exposure time and an aperture of 10 sq. arcsec. 
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